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Abstract. In this paper, we discuss the properties of CMDs, age, metallicity 
and radial velocities of eight massive LMC clusters using data taken from FORS2 
multiobject spectrograph at the 8.2-meter VLT/UT1. The strong near-infrared Ca 
II triplet (CaT) lines of RGB stars obtained from the high S/N spectra are used to 
determine the metallicity and radial velocity of cluster members. We report for the 
first time spectroscopically determined metallicity values for four clusters based on 
the mean [Fe/H] value of ~10 cluster members each. We found two concentrations 
in the distribution of ages of the target clusters. Six have ages between 0.8-2.2 Gyr 
and the other two, NGC 1754 and NGC 1786, are very old. The metallicity of the 
six intermediate age clusters, with a mean age of 1.5 Gyr, is —0.49 with a scatter of 
only 0.04. This tight distribution suggests that a close encounter between the LMC 
and SMC may have caused not only the restart of cluster formation in the LMC but 
the generation of the central bar. The metallicity for the two old clusters is similar to 
that of the other old, metal-poor LMC clusters. We find that the LMC cluster system 
exhibits disk-like rotation with no clusters appearing to have halo kinematics and 
there is no evidence of a metallicity gradient in the LMC, in contrast with the stellar 
population of the MW and M33, where the metallicity decreases as galactocentric 
distance increases. The LMC's stellar bar may be the factor responsible for the dilution 
of any kind of gradient in the LMC. 
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| Based on observations collected with the FORS2-VLT of the ESO within observing program 
078.D-0428(B). 
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1. Introduction 

The Large Magellanic cloud (LMC) is an important target in the effort to understand 
the stellar populations of galaxies, due to its close proximity, favorable viewing angle, 
ongoing star-formation activity (Harris & Zaritsky, 2009) and relatively low metallicity. 
The LMC cluster system is well known to show a puzzling age distribution, with a 
handful of old (~13 Gyr), metal-poor globular clusters and a number of intermediate- 
age (1-3 Gyr), relatively metal-rich populous clusters, with only a single cluster 'ESO 
121-SC03' in between (Mackey et al. 2006 and Mateo et al. 1986). The chemical 
enrichment of the LMC has been studied in some detail recently by Harris & Zaritsky 
(2009) and Carrera et al. (2008). The early history of the LMC was characterized by 
rapid enrichment, reaching [Fe/H]=— 1.2 by about 12 Gyr ago (Harris & Zaritsky, 2009). 
The metallicity then increased only slightly throughout the long quiescent epoch. When 
cluster formation resumed some 3 Gyr ago, the chemical abundances began a steady 
ascent to the present-day value of [Fe/H]=— 0.2 — 0.3 (Carrera et al. 2008; Hill et al. 
1995; and Rolleston, Trundle & Dufton 2002). 

These features provide clues to the history of star formation in the LMC. There 
was an initial epoch of star formation during which a significant portion of the LMCs 
stellar mass was formed. This era was followed by a long quiescent epoch, during 
which cluster formation was suppressed throughout the LMC. Then, about 3 Gyr ago, 
cluster formation resumed throughout the LMC, and has been ongoing since. The global 
resumption of star formation is strongly suggestive of a dramatic change at this point 
in the LMCs past, such as the merger of a gas-rich dwarf galaxy or a particularly 
dramatic tidal encounter, probably with the Small Magellanic cloud (SMC) given the 
similar resumption of star formation seen in that galaxys star formation history (SFH) 
(Harris & Zaritsky, 2009). Bekki et al. (2004) suggest that the recent burst of cluster 
formation is linked to the first very close encounter between the Clouds about 4 Gyr 
ago, which would have induced dramatic gas cloud collisions, allowing the LMC to begin 
a new epoch of cluster and star formation; strong tidal interactions between the Clouds 
have likely sustained the enhanced cluster formation. Bekki et al. (2004) also find that 
the close encounter between the LMC and SMC would have been sufficient to cause 
the formation of the LMC Bar around the time of the new epoch of cluster formation, 
giving rise to the similar SFHs seen in the cluster system and the bar (Cole et al. 2005, 
hereafter C05). In addition to enhancing star formation, tidal forces can result in the 
in-fall or outflow of material, thereby affecting the chemical evolution history (CEH) of 
the LMC and, at the same time, leaving behind a signature of the interaction. Thus, 
accurate knowledge of the ages and metallicities of LMC clusters is necessary to fully 
understand the formation and dynamical history of this galaxy. Here it is worthwhile 
to mention that recently Kallivayalil et al. (2006, 2009), using HST data for fields in 
the Magellanic Clouds (MCs) centered on background quasars, determined the proper 
motions of the LMC and SMC. They argued that the MCs may not be bound to each 
other or the Milky Way (MW) and may in fact be passing the MW for the first time, 



3 



although, to further verify and refine these results, more precise data is needed. 

Red giant branch (RGB) stars are among the brightest red stars in stellar systems 
that are older than ~ a Gyr. They appear in almost all galaxies. Therefore, red giants 
are an excellent tool for probing the parameters of old populations and the history of star 
formation in almost any galaxy. Globular clusters, with their single age and metallicity, 
are the ideal sites for calibrating the RGB parameters. Since Da Costa & Armandroff 
(1990) explored the possibility to use the color of the RGB as a metallicity indicator, 
many signicant advances have been made, both in the improvement of the astronomical 
instrumentation and the development of the corresponding theory. Kuchinski et al. 
(1995) and Minniti et al. (1995) demonstrated both empirically and theoretically that 
the slope of the RGB in a K vs. (J — K) color-magnitude diagram (CMD) is sensitive 
to the metallicity of the population. They investigated this correlation for metal-rich 
globular clusters and derived a linear relation between [Fe/H] and RGB slope. Ivanov 
et al. (2000) extended this relation to metallicities as low as [Fe/H] = —2.1. Ferraro et 
al. (2000) derived a similar relation. Ivanov & Borissova (2002) tabulated the relation 
in [M K , (J — Kg), [Fe/H]] space based on a large sample of 2MASS photometry of MW 
globular clusters. Tiede et al. (1997), using 4 Galactic open clusters, derived the slope 
of the RGB - metallicity relation for a younger population and Kyeong & Byun (2001) 
re-derived the calibration adding 10 new Galactic open clusters. And finally, Tiede et 
al. (2002) found two very clear linear trends, one for Galactic open clusters spanning 
—0.6 < [Fe/H] < +0.1 and another for Galactic globular clusters spanning —2.2 < 
[Fe/H] < -0.4. 

To derive the general relation between slope of the RGB, [Fe/H] and age, we have to 
extend the existing calibrations to the clusters in external galaxies, because there are no 
Galactic counterparts of the young metal-poor MC clusters. Accordingly, we observed 
more than a dozen clusters in the LMC and SMC in December, 2004 with the WFIRC 
attached to the 2.5-m du Pont telescope in the Las Campanas observatory. Thus, we 
have very good-quality CMDs, which allow us to determine the slope of the RGB of 
the clusters. The ages of the clusters are determined from good- quality HST CMDs 
with well-defined main-sequence turn-off points. The metallicities of most clusters in 
our RGB slope sample are based on the latest spectroscopic measurements available in 
the literature (GOG). There are eight clusters given in Table [H however, with no modern 
metallicity determination. There is an empirically developed, simply calibrated method 
available which can make an efficient estimate of metallicity ([Fe/H]) for individual 
RGB stars using the strength of the Ca II triplet (CaT) lines at 8498, 8542, and 8662 
A(Battaglia et al. 2008; Carrera et al. 2007; Cole et al. 2004, hereafter C04; C05; 
Grocholski 2006, hereafter G06; Parisi 2009, hereafter P09; Warren & Cole 2009 etc.). 
A powerful application of this method to the age-metallicity relation and abundance 
gradient of the star clusters of the LMC and SMC were already demonstrated by G06 
and P09, following up on the pioneering work by Olszewski et al. (1991, hereafter 
091). Therefore, to derive the metallicity, we obtained medium-resolution spectra at 
the near-infrared CaT using the FORS2 multiobject spectrograph at the 8.2-meter 
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VLT/UT1 with high signal-to-noise (20 < S/N < 100) ratio. The benefit of using 
CaT is that it is very efficient; not only are RGB stars near their brightest in the 
infrared, but the multiplexing capabilities of many moderate-resolution spectrograph's 
allow the observation of dozens of stars simultaneously, greatly increasing the probability 
of identifying cluster members (P09). In addition, the metallicity calibration is very 
straightforward and well calibrated and independent of luminosity, distance, reddening 
and age. Therefore, using FORS2 on the VLT, we were able to identify 77 member 
stars in eight massive LMC clusters and determine accurate mean cluster velocities and 
metallicities. 

In this paper, we discuss the metallicity distribution, their relationship and the 
kinematics of LMC clusters, extending the work done by G06. We will further discuss 
the relation between RGB slope, metallicity and age in an upcoming paper (Kurtev et 
al. 2010). In §2, §3, §4, and §5, we discuss the observations, techniques required to 
determine radial velocity & metallicity and the cluster membership criteria. In §6 & §7, 
we present the basic cluster parameters based on our photometry and spectroscopy and 
analyze the results, followed by conclusions in §8. 

2. Observations and data reduction 

2.1. Photometry 

Short exposure images of eight LMC clusters were observed with FORS2, VLT, through 
the filters V and R under excellent atmospheric conditions (seeing < 0.7 arcsec) on 5 
October 2006 and 17 November 2006. FORS2 was used in normal resolution with the 
pixel size of 0.25"/pixel and field of view of 6.8' x 6.8', which allowed us to sample, in the 
outermost field regions, a good portion of the background LMC field. The coordinates 
of the clusters as well as the total exposure times are listed in Table [TJ 

De-biasing and flat-fielding were performed with standard IRAF packages, using 
the calibration frames (bias and sky flat fields) obtained during the same nights 
as the target observations. Photometric reduction was carried out with the 
DAOPHOT II/ALLFRAME package (Stetson 1987, 1994). Stellar PSF models for each 
frame were obtained by means of a large set of uncrowded and unsaturated stars. The 
mean photometric errors are 0.003-0.009 mag for (V,R) < 20 mag and between 0.01 
mag and 0.03 mag for the fainter magnitudes. The nights were photometric and the 
calibration to the Johnson standard photometric system was performed by comparison 
of the stars in the Landolt (1992) field PG2213-006. The maximum zero-point errors 
are calculated as ± 0.01-0.02 magnitudes. 

2.2. Spectroscopy 

The spectroscopic observations were carried out with FORS2 in visitor mode at the 
Antu (VLT-UT1) 8.2 m telescope at ESOs Paranal Observatory during the night of 16- 
17 November 2006. We used the FORS2 spectrograph in mask exchange unit (MXU) 
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mode, with the 1028z+29 grism and OG590+32 order blocking filter. The MXU slit 
mask configuration allows the placement of more slits (in this case 37) on the sky 
than the 19 movable slits provided in Multi Object Spectrograph mode. Targets were 
selected so as to maximize the number of likely cluster members observed; typically 
20 stars inside our estimated cluster radius (see §5) were observed, with an additional 
15 stars outside of this radius that appeared to be LMC field red giants based on our 
infrared photometry with the du Pont telescope. 

We used slits that were 1" wide and 6.25" long. In all fields, we obtained 3 to 4 
separate 600s exposures, taken under seeing less than 1". They were average combined to 
form the final spectra. Pixels were binned 2x2, yielding a plate scale of 0.25" pixel -1 . 
We used the upper (master CCD) of FORS2 for observations of the target clusters 
which has a readout noise of 2.9 e - and inverse gain of 0.7 e - ADU -1 . The spectra 
have a dispersion of ~ 0.85 A/pixel (which corresponds to a resolution of 2-3 A) with a 
characteristic rms scatter of 0.06 A, and cover a range of ~ 1600 A in the region of the 
CaT (8498, 8542, and 8662 A). For the final spectra, S/N is typically 40-60 pixel" 1 with 
some stars as high as ~ 100 pixel -1 and, in some few cases, as low as ~ 20 pixel -1 . 
For the present study we use spectra having S/N > 35 pixel -1 . Calibration exposures, 
bias frames and flat-fields were also taken as a part of the standard calibration plan. 

We used the IRAF task ccdproc to fit and subtract the over-scan region, trim the 
images, and flat-field each image. Flats were corrected for the lamp shape by the task 
response. We used the task apall to define the sky background and extract the stellar 
spectra onto one dimension. The sky level was defined by performing a linear fit across 
the dispersion direction to sky windows on each side of the star. When the star fell near 
the top or bottom of the slitlet; the sky regions were chosen interactively. The tasks 
identify and dispcor were used to calculate and apply the dispersion solution for each 
spectrum. Finally, the spectra were continuum normalized by fitting a polynomial to 
the stellar continuum. Examples of the final spectra in the CaT region can be seen in 
Figure [TJ Clearly, the three Ca II lines dominate the red giants spectra in this region. 
Also note the stronger Ca lines in the bottom spectrum, which, given the very similar 
T e ff and log g values (as indicated by their virtually identical V — Vhb values, i.e. the 
mag difference between the star and the horizontal branch), graphically illustrates the 
higher metallicity of the bottom spectrum and the power of the CaT technique. 

3. Radial Velocity 

Radial velocity measurements for our target stars has been done to help us in defining 
the cluster membership as well as to study the kinematics of the LMC. Since a clusters 
velocity dispersion is expected to be relatively small compared to the surrounding field 
and its mean velocity quite possibly distinct from the field, radial velocities are an 
excellent tool for determining cluster membership. 

Radial velocities for all the target stars were determined through cross-correlation 
with eight template spectra from NGC 1851 (Carerra et al. 2007) using the IRAF task 
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fxcor (Tonry & Davis 1979). These template stars were observed as a part of the CaT 
calibration program; thus, this sample offers a good match to the spectral types of our 
target stars. In addition, these observations were made with a telescope and instrument 
setup that is almost identical to those used by us. In addition to calculating relative 
radial velocities, fxcor uses information about the observatory location, date and the 
time of the observations (once the ESO header has been appropriately reformatted) 
to correct the derived velocities to the heliocentric reference frame. For a stars final 
heliocentric radial velocity, we adopt the mean of all cross-correlation results. We find 
good agreement among the template-derived velocities, with a typical standard deviation 
of 8 km s" 1 . Having used the same set up as G06, we follow the same approach to 
calculate the error due to imprecise alignment of the slit center and the stellar centroid. 
G06 estimate this latter error as ±4.2 km s _1 . This, together with the the standard 
deviation of 8 km s _1 gives a typical error of 9 km s _1 in measuring the radial velocity 
in the present study. 

4. Equivalent widths and Abundances 

The abundances of the clusters in the present sample have been derived from the 
equivalent width (EW) of strong CaT lines in individual RGB stars based on the 
empirically developed method which has been explored already in detail previously 
(see for eg. C04, C05, G06, P06). The CaT lines can be contaminated by weak neutral 
metal lines, and the continuum may be affected by weak molecular bands. Both effects 
make it impossible to measure the 'true' EWs of the CaT lines in isolation. Armandroff 
& Zinn (1988) defined two pseudocontinum windows in line free regions on either side 
of each feature. For each of the Ca II lines, the line connecting the average number of 
the counts in each of these two windows was defined to be the pseudocontinum. The 
pseudo-equivalent width was then computed by taking the integral over the feature 
boundaries of the difference between the pseudocontinum and the spectrum. Further, 
C04 mentioned that the spectra need not be flux-calibrated because the continuum 
slopes of red giants in this wavelength range are generally close to flat, so a simple 
normalization by an arbitrary polynomial suffices for measurement of the pseudo- 
equivalent width. In the present study, we calculate the EW by direct pixel-to-pixel 
integration using the interactive splot program of IRAF. The EWs are measured by 
visually marking two continuum points at the line band pass given by Armandroff & 
Zinn (1988). The continuum level, which is a linear function between these two marked 
points, is effectively fixed at 1 by our normalization procedure. The integral is taken 
between these two points. Since our spectra are of high enough quality that all three 
CaT lines are well measured, we define our metallicity index Y,W, as the sum of the 
EWs of three CaT lines (see for eg. C04, G06, and P09). 

Although theoretical and empirical studies have shown that effective temperature, 
surface gravity, and [Fe/H] all play a role in CaT line strengths (e.g., Jorgensen et al. 
1992; Cenarro et al. 2002), Armandroff & Da Costa (1991) showed that there is a 
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linear relationship between a stars absolute magnitude and TW for red giants of a given 
metallicity. Following previous authors using the same CaT technique (G06 & P09), 
we define a reduced equivalent width, W, to remove the effects of surface gravity and 
temperature on TW via its luminosity dependence: 
W = ?W + P(V - V HB ), 

The introduction of the difference between the visual magnitude of the star (V) 
and the clusters horizontal branch/red clump (V H b) in calculating W, also removes 
any dependence on cluster distance and interstellar reddening. The value of (5 has been 
investigated in detail by previous authors (C04, G06, P09) and we have adopted the 
value obtained by C04: (5 = 0.73 ± 0.04. We chose C04s value of (3 because they used 
both globular and open clusters, covering a broad range of ages, to derive it and their 
instrumental setup is very similar to ours. We note that, as discussed in detail in G06, 
the use of C04s value for (3, combined with defining the brightness of our target RGB 
stars relative to the HB of their parent cluster, incorporate any age effects into the CaT 
calibration, thus it is not necessary to make any corrections for the ages of our clusters. 

As shown by Rutledge et al. (1997), there is a linear relationship between a clusters 
reduced EW and its metallicity on the Carretta & Gratton (1997) abundance scale for 
globular clusters of the MW. C04 extended this calibration to a wider range of ages 
(2.5 Gyr < age < 13 Gyr) and metallicities (—2.0 < [Fe/H] < —0.2) by combining the 
metallicity scales of Carretta & Gratton (1997) and Friel et al. (2002) for globular and 
open clusters, respectively. Although many of our clusters are younger than 2.5 Gyr, 
we adopted the C04 relationship, 

[Fe/H] = (-2.966 ± 0.032) + (0.362 ± 0.014)W, 

to derive the metallicities of our entire cluster sample. Carrera et al. (2007) 
investigated the behavior of CaT in the age range of 0.25 Gyr < age < 13 Gyr and 
the metallicity range —2.2 < [Fe/H] < +0.47. In spite of the extended age range 
of their calibration, we decided to use C04s calibration for three reasons. First, the 
calibration of Carrera et al. (2007) uses the absolute magnitude of a target star (rather 
than the brightness relative to the HB, as we have adopted), thus requiring accurate 
photometric calibration and distance to the star. Second, using (V — Vhb) instead of 
absolute magnitude removes the need to make any assumptions about the foreground 
reddening toward a cluster. Finally, Carrera et al. (2007) performed estimates of the 
differences in W as a function of age using the Jorgensen et al. (1992) models and BaSTI 
stellar evolution models (Pietrinferni et al. 2004). From their data, they confirmed that 
the influence of age is small (for clusters having age 10.5 and 0.6 Gyr, the difference in 
metallicity is ~ 0.25 — 0.15 dex), even for ages < 1 Gyr. Based on these results, we are 
confident that the C04 calibration can be accurately extended to younger ages (see also 
P09). 

In addition, Battaglia et al. (2007) have shown that [Fe/H] values derived from 
CaT agree with those from high resolution spectroscopy and detailed model atmosphere 
analysis to within 0.1 - 0.2 dex over the range -2.5 < [Fe/H] < -0.5, based on FLAMES 
data for a large sample of stars in the Sculptor and Fornax dwarf spheroidal galaxies. 
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It is well known that such stars have distinct [Ca/Fe] ratios from those in the Galactic 
globular cluster calibrators (e.g., Shetrone et al. 2001; Geisler et al. 2007) and also 
possess a larger range of stellar ages, extending to much smaller values. Nevertheless, 
CaT yields metallicities that are very close to real Fe abundances, reinforcing our 
confidence in the technique. 

With the assumption that each star in a given cluster has the same metallicity and 
thus the same W value, the error on W, aw', depends on the standard deviation of 
the individual star's W values about the mean cluster value as well as the number of 
stars, N, used in each cluster (see also Warren & Cole 2009) i.e. aw' = —r — , this 



error is further propagated to define errors in calculating the [Fe/H]. Our estimate of 
the systematic error per star ranges from 0.08 to 0.14 dex, with a mean of 0.11 dex. 

5. Cluster Membership 

Cluster membership criteria are very crucial in defining the parameters for the clusters. 
We use a combination of several criteria to isolate cluster members from field stars as 
defined by G06. First, all observed targets were selected on the base of the statistically 
decontaminated near-IR CMDs (Kurtev et al., 2010). Second, the cluster centers and 
radii are chosen by eye, based primarily on the photometric catalog (VR). As an 
example, Figure [2] (left panel) shows (x, y)-positions for all stars photometered in the 
NGC 1754 and NGC 1900 fields respectively, with large dots denoting our target stars 
and the large open circle representing the adopted cluster radius; stars outside this circle 
are considered non-members due to their large distance from the cluster center. We note 
that stars outside the cluster radius were observed so as to define parameters for the LMC 
field, which aids in isolating cluster members. Next, radial velocity versus radial distance 
is plotted in Figure [2] (right panel). Stars moving with the velocity of the clusters are 
easily identified due to their smaller velocity dispersion. Our velocity cut, denoted by the 
horizontal lines, has been chosen to include the expected observed velocity dispersion in 
each cluster. To determine this, we have adopted an intrinsic cluster velocity dispersion 
of 5 km s -1 and added this in quadrature with our adopted radial velocity error, 9 km 

resulting in a dispersion of ~ 10.3 km s -1 . Thus, we adopted a width of ±11 km 
for our radial velocity cut. We select the cluster velocity to maximize the number of 
stars lying within this range. The cluster radius (Figure [2] (right panel), vertical line) is 
also marked for reference. Finally, Figure [3] shows metallicity as a function of distance 
for the stars in NGC 1754 and NGC 1900 respectively, with horizontal lines representing 
the [Fe/H] cut of ±0.2 dex that has been applied to these data. We adopt random error 
of ±0.15 dex in measuring each individual [Fe/H], this along with a typical systematic 
error of ±0.11 dex gives total error in measurement ±0.18 dex. Thus, we have rounded 
this figure to ±0.2 dex. Radial velocity members which are in this range of metallicity 
are assumed to be probable cluster member. Star symbols denote stars that have made 
all of these cuts and are therefore considered to be cluster members. 

In Figure H] we present the traditional SW versus V — Vhb plot for cluster members 
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of NGC 1754 and 1900 with the lines representing the mean metallicity of —1.48 and 
—0.55, respectively, with a standard error of the mean =0.03. The CMD in Figure [5] 
shows all stars photometered along with CaT targets (big dots) in the NGC 1754 and 
NGC 1900 field respectively. Confirmed member RGB stars are also shown by star 
symbols. 

In Table [2], for all RGB stars of the observed LMC clusters, we list the following 
information: stellar identification number, RA, Dec, heliocentric radial velocity and its 
associated error, V — Vhb, S/N, EW, W along with errors, and the value of [Fe/H] with 
associated error. 

6. Cluster Properties 

6.1. Color- Magnitude Diagrams 

Figure [6] shows the CMDs of the eight LMC clusters. All the stars measured with good 
precision (error in mag <0.1) in the whole FORS2 frames are plotted here: therefore the 
CMDs include both cluster and field stars. To check the zero point of our calibrations 
we over-plotted the fiducial lines of RGB stars, field RR Lyrae and Red Clump (RC) 
stars as defined from MACHO database (Alcock et al. 2000, see their Table 1.). 

The clusters from our sample are located in areas densely populated by field stars. 
In order to obtain a cleaner cluster CMD, we proceeded in the following way. For each 
cluster, an annulus around the center was selected with the criterion of being large 
enough to contain most cluster stars, but small enough to minimize the background 
contribution. The small circular region containing the very center of some of the clusters 
was excluded because these areas were too crowded, and therefore the photometry was 
rather poor. The photometry around very bright stars projected close to the objects 
NGC 1786, SL826, and NGC 1754 is also removed. A CMD typical of LMC field stars 
was extracted from a second region, well outside a given radius around the cluster 
center, selected with the criterion of containing a number of stars large enough to insure 
good statistics, but far enough from the cluster not to contain a significant fraction 
of cluster stars. The cluster and reference fields cover the same area. The full range 
of magnitude and colors of a given CMD was divided into grids and we computed the 
expected number-density of field stars in each cell based on the number of comparison 
field stars. After that the expected number of field stars is randomly subtracted from 
each cell. The statistically de-contaminated CMDs are plotted in Figure [6] as dark dots. 

To check our de- contamination procedure we calculated the stellar surface-density 
distributions (number of stars per arcmin 2 ) (Bonatto & Bica, 2006). They are computed 
for a mesh size of 0.2 x 0.2 arcmin 2 , centered on the apparent cluster center. Figure[7J 
left panel shows observed photometry of NGC 1795, while on the right-hand panel only 
the stars selected by means of the color-magnitude filter described above are plotted. 
As can be seen, most of the field stars have been removed, although a small number of 
field stars still fall in the cluster main sequence locus region. 
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The decontaminated CMDs outline the basic features of the clusters. Most of them 
show typical diagrams of intermediate age clusters. As can be seen in Figure |8] we can 
not reach the Turn off/Termination Point (TO) for NGC1795; NGC1754; NGC1786; 
Hodge 2 because of the crowding and short exposure times. For the rest of the clusters 
from our sample: NGC1900; SL509; SL817; and SL862 the main sequence can be 
followed up to 4 magnitudes bellow the TO. The RC stars are located between « 19- 
19.3 mag and the RGB extends between 2 and 4 magnitudes brightward of the RC. 
NGC1754, and to a lesser extent NGC1786, show well populated horizontal branches 
(HB). 

6.2. Age determination 

The goal of this project is to obtain homogeneous [Fe/H] values and ages to be used 
for calibration of the near-infrared _fT s -band RGB slope-metallicity-age relation. The 
determination of the metallicities is described in §4. The most reliable technique for age 
determination is to examine the luminosity of the main sequence turnoff and compare it 
to the predictions of theoretical isochrones. Unfortunately, it is not possible for half of 
our objects, since the photometry is just not deep enough and we could not reach the TO 
and main sequence region. We choose to compare the rest of our CMD with the Padova 
set of theoretical isochrones (Marigo et al. 2008), because they are based on models 
with a moderate amount of convective overshooting, asymptotic giant branch evolution 
and circumstellar dust around the coolest evolved stars. For the clusters NGC 1795, 
NGC 1754 and NGC 1786 much deeper, better HST data are available in the literature 
and we adopted the age determination from the corresponding papers. Since in these 
papers different sets of theoretical isochrones are used, we plotted the corresponding 
isochrone from the Padova set to check the coincidence. To determine the best models 
we applied a statistical tool to compare the observed and theoretical main sequence 
and/or red clump. We compared the differences along the observed and theoretical 
fiducial lines of each cluster using a £ 2 statistics. 

NGC 1795: Because of the crowding, incompleteness and short exposure times of 
our images the CMD of NGC 1795 does not reach the TO. The RGB and RC stars 
however are clearly visible and the corresponding branches are well populated. To 
determine the RC brightness we adopted the method described in Pietrzyhski et al. 
(2003): we selected the stars with V-band magnitudes and V — R colors in the range 
around the RC stars and fit the derived l/-band magnitude histogram with a complex 
function involving a Gaussian function (representing the distribution of RC stars) 
superimposed on a second-order polynomial component (for the stellar background; 
see, e.g., Stanek & Garnavich 1998; Alves et al. 2002; Pietrzynski & Gieren 2002). 
Such a described procedure finds the RC at V A =19.13±0.07. The latest photometric 
investigation based on the archival HST images of the cluster is given in Milone et al. 
(2009). The observed deep (V, V — I) CMD is compared with a set of BaSTI isochrones 
and gives the following parameters: (M — m) =18.45; E(B — V)=0.10, z=0.008 and age 
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1.3 Gyr. Since we have accurately determined metallicity [Fe/H]=-0.47 we adopted their 
distance, reddening and age values and plot the corresponding isochrone from Padova 
set which fit well the luminosity and color of the RC region (Figure El upper left panel). 

NGC1754: NGC 1754 is well known old LMC cluster (Olsen et al. 1998). The main 
features visible on our CMD are well populated RGB, long blue HB, typical for metal 
poor clusters and a sub-giant branch (SGB) down to V=22 mag. The TO point and 
RC regions are not visible. In our photometry, the RC is determined at 1/=19.25±0.17. 
The level of the HB calculated as lower envelope is found at V A =19.55±0.11. Olsen et al. 
(1998) derived (M — m) =19.00 (based on the MW globular cluster stars comparison) 
and (M - m) =18.62 (based on the level of the HB stars); E(B - V)=0.10 and age 
15.6 Gyr. We adopted their average distance value of (M — m) =18.66 and plotted the 
Padova's isochrone of 15.6 Gyr, but the deviation from the observed CMD is very large. 
Taking into account our metallicity [Fe/H]=— 1.48 the best fit of the Padova isochrones 
(Marigo et al. 2008) favorite an age of 14 Gyr. To check this age determination we 
retrieve the HST V, I photometry of Olsen et al. (1998). As can be seen in Figure El 
(right corner of NGC 1754 plot) an age of 14 Gyr fits well the TO region of the Olsen 
et al. (1998) CMD. Obviously, this is slightly older than the well-determined age of the 
Universe from WMAP and therefore an overestimate. 

NGC 1786: This object is another example of an old, metal poor cluster in the LMC 
(Brocato el al. 1996). In our CMD the RGB is well defined, the HB is populated on 
the blue side of the instability trip at y=19.44±0.08 (lower envelop). The fundamental 
properties of NGC 1786 are taken from Table 7 of Grocholski et al. (2007) and are based 
on the mean apparent l^-band magnitude of RR Lyrae stars (Walker 1985; Walker & 
Mack 1988; Walker 1989, 1990, 1992ab, 1993). With (M — m) =18.62; E(B - V)=0.06 
and [Fe/H]=-1.77, we plotted Brocato el al. (1996) isochrone of 12.3 Gyr. As can be 
seen in Figure [8] the isochrone for the adopted age fits well our CMD. 

Hodge 2: There is no photometric information for this cluster in the literature. Our 
CMD does not clearly reach the TO point; however, the RGB is well defined and the 
RC can be located at Vrc= 19.12±0.17. To determine the age of the cluster we need 
reddening and distance. We calculate the distance using the Population II Cepheid star 
SV*HV 2439 (OGLE LMC-SC7 21841), a confirmed cluster member (Hodge & Wright 
1963; Udalski et al. 1999). According to OGLE III database (Soszynski et al. 2008) 
the apparent mean V magnitude of the star is V=15. 58, with a period of pulsation of 
P=4.81308 days. The reddening E(B - V)=0.U is adopted from the OGLE database 
as a mean value of the SC7 field. The distance modulus (M — m) =18.62 is calculated 
as a difference between reddening corrected mean V magnitude of SV*HV 2439 and the 
absolute My magnitude of Cepheid stars given in Benedict et al. (2002). The cluster 
has a metallicity value [Fe/H]=— 0.49 ± 0.11 and the best fit isochrone favors an age of 
1.6 Gyr. 

NGC 1900: The CMD presented here is the only photometric investigation up to 
the moment and shows typical features of an intermediate age stellar cluster: long and 
well populated main sequence reaching about 4 magnitudes below the turnoff point and 
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some evolved RGB stars. The distance of (M — m) =18.47±0.07 is calculated from 
the mean RC K$ brightness derived from our high quality J, K$ infrared photometry 
(Kurtev et al. 2010). We used the absolute RC magnitude Mk = —1.644, calibrated in 
Alves et al. (2002). To calculate the distance a correction for population effects has to 
be applied (Gullieuszik et al. 2007), because of the differences in the stellar content of 
LMC and Galactic RC stars on which the Alves (2002) calibration is based. Salaris & 
Girardi (2002) estimated the population effects on the RC absolute magnitude in the 
K band of AM§ C = -0.03. Using [Fe/H]=-0.55 ± 0.09 obtained in §5, we derive an 
age of 800 Myr. 

SL 509: The only photometric investigation up to the moment is reported by Bica et 
al. (1998) and Piatti et al. (2003) and is based on the Washington photometric system. 
The derived age in Piatti et al. (2003) is 1.2 Gyr and the metallicity corrected for age 
effects is [Fe/H] =— 0.65. They also reported a secondary clump 0.45 mag fainter and 
bluer than the main RC, which also can be noted in our CMD. In Figure [9] (left panel) we 
zoomed the RC and TO regions of the CMD, marking with boxes the suspected clumps. 
The luminosity function of the RC region is shown in the middle panel of Figure [9] and 
finally the last panel of the Figure [9] represents the positions of the stars with respect to 
the cluster center. As can be seen there is a strong hint of two populations, as pointed 
out by Piatti et al. (2003), but the sample is too small with only 66 stars and more 
data are necessary to confirm the reality of these populations. If confirmed, according 
to Girardi et al. (2009) the secondary clump could be a result of He-ignition in stars 
just massive enough to avoid e~-degeneracy in their H-exhausted cores. The TO can 
be located at Vto= 19.40±0.11. Using the same method as for NGC 1900 we derived 
(M-m) =18.48±0.05. The E(B-V)=0.03 is taken from the OGLE database. The best 
isochrone fit gives 1.2 Gyr, in excellent agreement with Piatti et al. (2003) result. 

SL817: The cluster has only Washington photometry (Bica et al. 1998; Piatti et 
al. 2003) and shows morphology similar to that of SL509. Using the same technique 
as for SL509 we derived the TO at Vto= 20.20±0.11, the RC is calculated at Vrc= 
18.96±0.17. For the distance, reddening and age we determined (M — m)o=18.38±0.07; 
E(B-V)=0.12 and 2.2 Gyr, respectively. 

SL862: And finally for SL862 we have TO at V TO = 20.10±0.11; the RC at V RC = 
19.07±0.09; (M - m) =18.45±0.07; E(B-V)=0.10 and an age of 1.7 Gyr. Piatti et al. 
(2003) give an age of 1.7 Gyr and the metallicity of —0.75 from Washington photometry, 
in agreement with our values. 

We estimated the errors in the age determination as an upper limit of the various 
factors: errors in the photometry, errors in the transformation to the standard system, 
errors due to incompleteness, field decontamination, metallicity determination and errors 
of the fit to the theoretical isochrones. For the clusters without good main sequence 
photometry, ages can be estimated to within ± 1.4-1.7 Gyr, while for the younger 
clusters with prominent main sequences we estimate an error of ±0.5 Gyr. In Tabled 
we have listed the obtained value of age for all the target clusters. The radial velocity 
and metallicity along with standard deviation is also given in Table [3j which has been 
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obtained by taking the mean of the individual values of cluster members. 
7. Analysis 

7.1. Comparison to previous metallicity determination and sample enlargement 

We compared the values of the metallicities obtained in the present work with previously 
published values. Since our method relies on high S/N spectroscopy of a large number 
of stars per cluster, the individual and mean errors are in general smaller as compared 
to the previously determined values (for eg. 091 have used in general 2 stars per cluster, 
while we have on average 10 identified cluster members). Four of our target clusters have 
metallicity values given in 091 using CaT spectroscopy of some individual giants. Five 
of our target clusters also have photometrically determined values of metallicity (see 
Table [3]). As indicated by G06, although 091 used CaT lines as a proxy for measuring 
Fe abundance directly and is thus similar to our work, it is inappropriate to directly 
compare the 091 values to our work due to differences in the details of the metallicity 
determination (see for detail G06). C05 find that one can estimate the abundance of 
091 clusters on the metallicity system we have used via the following conversion: 
[Fe/H] ~ - 0.212 + 0.498 [Fe/H] 09 i - 0.1 28 [Fe/H] ^ 91 

In Table El column (6), (7), (9) gives the [Fe/H] values obtained in the present 
study, values given by 091 (spectroscopically) and the photometrically determined 
values, respectively. In Column (8), we give the transformed values of [Fe/H] given by 
091 to our system using the above relation. In Fig. [TDl we have plotted the difference 
£ [Fe/H] presen< - [Fe/H] Kterature ' as a function of age for 091 clusters (filled circles: original 
091 values, open circles: transformed 091 values) and photometrically observed clusters 
(open triangles). 

The difference of metallicity between the present value and the transformed value 
of 091 are of the order of 0.15 — 0.23 dex. 091 give their [Fe/H] errors for individual 
stars as 0.2 dex, therefore, deviations between these data sets as large as ±0.2 are not 
unexpected, suggesting that these results are in relative agreement with no offset. We 
note, however, that even with the use of the above equation, it is very difficult to directly 
compare the derived cluster abundances because of the difference in calibration strategy. 
Comparison of our [Fe/H] values with the photometrically determined values for each 
cluster is given below: 

NGC 1754: This cluster has been studied by Olsen et al. (1998) using HST 
CMDs quoting a metallicity of —1.42 ± 0.15. Olsen et al. (1998) commented that the 
abundances they obtained are on average higher than previously measured spectroscopic 
abundances. They have used the method of Sarajedini (1994) to determine abundances 
from the measurement of the height of the RGB above the level of HB. 

NGC 1786: This cluster has been studied by Brocato et al. (1996) using VI data 
from the 3.6 ESO NTT. They have determined a lower value of —2.1 but with a large 
error (0.3 dex) using the Sarajedini (1994) method. Recently this cluster was also 
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studied by Mucciarelli et al. (2009) using high-resolution spectroscopy of giant stars. 
Their value of [Fe/H]= — 1.75±0.01 is virtually identical to our metallicity determination 
([Fe/H] =—1.77 ± 0.08), lending confidence to our results. 

SL509, SL817 & SL862: These clusters have been studied by Piatti et al. (2003) 
using Washington photometry. Our [Fe/H] values for SL509 and SL817 (i.e. — 0.54±0.09, 
—0.41 ± 0.05) are comparable to their reported values of —0.65 ± 0.2 and —0.35 ± 0.2, 
respectively. For SL862, they gave [Fe/H] = —0.75 ±0.2, which is lower than the present 
estimation of -0.47 ± 0.07. 

For all eight clusters, we have used the very reliable CaT method as the [Fe/H] 
indicator, with significantly more cluster members and less scatter in individual [Fe/H] 
value compared to previous studies. 

It is important to have significant numbers of data points to make a definitive 
analysis of the age-metallicity relation of LMC clusters. Over the past few decades, 
nearly a hundred LMC star clusters have had both accurate age estimated from isochrone 
fitting to their CMDs (Sagar & Pandey 1989; Balbinot et al. 2009; and references 
therein), and a metallicity measurement from either spectroscopy and/or isochrone 
fitting (G06). Recently Harris & Zaritsky (2009) have compiled the age and metallicity 
of 85 clusters from their literature search. Out of these 85, 5 are in common with our 
sample and 15 are in common with G06. We have excluded our 5 and 15 G06 clusters 
from this sample and generate a new sample of 101 clusters [(85 -5 -15) + 8 + 28 = 101]. 
The present work combined with G06 makes a homogeneous dataset of 36 clusters using 
the same method for metallicity determination. In the case of G06 clusters, out of 28, 
15 have ages from the sample of Harris & Zaritsky (2009). For the other 13 clusters, we 
converted the given SWB type (Searle, Wilkinson, & Bagnuolo 1980, hereafter SWB) 
into age from the conversion given in Bica et al. (1996). These 13 clusters have either 
SWB type V or VI. Based on published data, we have taken the mean age of SWB 
type V and VI as 1.5 and 2.5 Gyr respectively. Here it is important to mention that 
the Harris & Zaritsky (2009) compilation is heterogeneous, using a variety of different 
age and metallicity scales, so in fact it is a bit dangerous lumping them all together, 
although this effect is minimal for the present study. 

1.2. Age /Metallicity Distribution and their Relationship 

The cluster metallicity distribution (hereafter MD) is an important diagnostic of the 
global chemical evolution of a galaxy and is useful for an overall comparison of the 
cluster systems of different galaxies (P09). 

In Figure [11] we have shown first the distribution of the ages of all 101 clusters. 
The black, grey and open regions represent the data from the present study, from G06 
and from Harris & Zaritski (2009), respectively. Broadly we can classify the population 
in three age bins: < 1, 1 — 3 and > 12 Gyr, ie. youngest, intermediate age and old age 
populations. In Table HI we have presented the statistics for the value of [Fe/H] and 
age in these bins along with the number of clusters using the CaT method or clusters 
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from the literature survey. We can easily see the change in the value of [Fe/H] with age. 
Globally, for clusters having mean age ~ 0.2 Gyr, the [Fe/H] value is —0.35, for age ~ 2 
Gyr, the [Fe/H] value is —0.55 and for age ~ 14 Gyr, the [Fe/H] value is —1.62. This 
change is more prominent in 2 nd and 3 rd bin of clusters having used the CaT method 
(we have only one cluster in 1 st bin). We found that the younger clusters are generally 
more enriched than older clusters with the youngest bin [Fe/H] value comparable to 
the present-day LMC value of '—0.31 ± 0.04' given by Rolleston, Trundle & Dufton 
(2002) and '—0.27 ± 0.06' given by Hill et al. (1995). The MD is shown in Figure [El 
where the color code is the same as Figure [Til We can easily see the strong peak in the 
distribution of clusters having [Fe/H] = —0.5 to —0.4 and rest of them are distributed 
approximately uniformly with metallicity. The metallicity of the six young clusters 
obtained in the present study, having a mean age of 1.5 Gyr, is —0.49 with a scatter of 
0.04, in excellent agreement with the previous finding of G06, which shows a very tight 
distribution with mean metallicity of —0.49 (cr = 0.09) for the younger clusters, which 
formed after the age gap ended some 3 Gyr ago. The narrow metallicity distribution 
of LMC young-intermediate age clusters has also been observed previously, based on 
high resolution spectroscopy (i.e. NGC 1866: [Fe/H] = -0.5 ±0.1 (Hill et al. 2000) and 
mean [Fe/H] = -0.38±0.09 for 4 LMC globular clusters (Mucciarelli et al. 2008)). This 
agrees well with the theoretical work of Bekki et al. (2004), which indicates that a close 
encounter between the LMC and SMC caused not only the restart of cluster formation in 
the LMC but the generation of the central bar (mean [Fe/H] = — 0.37±0.15; C05). We 
can easily see the two separate distributions in the age of the clusters (Figure [UJ and 
older clusters are showing a much larger spread in [Fe/H] value as compared to younger 
ones (see also Figure [T2l) . The sample size is near the limit of drawing statistically 
significant results of this kind. It is not clear why the SMC cluster age and metallicity 
distribution are distinct from that of the LMC. It is of great interest to enlarge the 
present homogeneous sample in order to ascertain the LMC cluster MD. Note that the 
CaT data shows a much tighter metallicity distribution than that of Harris & Zaritsky 
(2009) or 091. 

The age-metallicity relation (AMR) is a record of the progressive chemical 
enrichment of the star-forming local interstellar medium during the evolution of the 
galaxy, and so provides useful clues about the star formation and chemical evolution 
history of the local environment (Carraro et al. 1998). All existing AMR studies 
suffer from uncertainties due to either inhomogeneous datasets or photometry based 
abundances. In Figure [131 we make an updated AMR using homogeneous set of data 
using the CaT method for [Fe/H] measurement (36 clusters, 8 from the present work 
(star symbols) & 28 from G06 (closed circles)). We have also plotted for comparison, 
the data from the sample of Harris & Zaritsky (2009) (open circles) which are not in 
common with either the present work or the work done by G06. We can easily see there 
is less scatter in [Fe/H] values for the clusters using the CaT method. 

The essential features of the AMR have not changed (see also Harris & Zaritsky 
2009): there is still a long age gap, which is also a metallicity gap, especially if one only 
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looks at the CaT data. Prior to the age gap, the LMC formed a number of low-metallicity 
populous clusters, analogous to the MWs globular cluster population. Following the age 
gap, cluster formation resumed in the LMC, and proceeded roughly continuously to the 
present day. The clusters formed following the age gap have less dispersion in metallicity 
and are signicantly more enriched than the ancient cluster population, and there is some 
evidence for further enrichment with time for ages < 1 Gyr. 

C05 have obtained stellar metallicities for almost 400 stars in the bar of the LMC, 
also using CaT lines. Their obtained AMR shows a similar behavior to that of the 
clusters for the oldest population. However, while for the clusters the metallicity has 
increased over the last 2 Gyr, this has not happened in the bar. Recently, Carrera et al. 
(2008) have examined the AMR for field stellar populations, based on CaT spectroscopy 
of individual red giants in four LMC fields. They find that the disk AMR is similar to 
that of the LMC star clusters: there was rapid initial enrichment to [Fe/H]~ —0.8 by 
about 10 Gyr ago. Following this, the metallicity increased only slowly until 3 Gyr 
ago, when it began a more rapid increase to its present-day value of [Fe/H]~ —0.2. 
However, the lack of objects with ages between 3 and 10 Gyr is not observed in the field 
population. Their reconstructed SFH (for their inner disk fields) indicates there was a 
relatively quiescent epoch between 5 and 10 Gyr ago but some star formation occurred. 
These conclusions differ in detail from those derived from the AMR of the central region 
of the LMC (C05), but the observed AMRs are entirely consistent. 

Carrera et al. (2008) also showed that the bar AMR differs from that of the disk 
in the last 3 Gyr: while in the disk the metallicity has increased during this time, in 
the bar it has remained approximately constant. This would be in agreement with the 
prediction that a typical bar instability pushes the gas of the disk towards the center 
(Sellwood & Wilkinson 1993), so the infall gas is expected to be previously pre-enriched 
by the disk chemical evolution. This would be also the case in the scenario suggested 
by Bekki & Chiba (2005), in which the bar would have formed from disk material as 
a consequence of tidal interactions between the LMC, the SMC and the MW about 5 
Gyr ago, the epoch when the bar AMR differed from that of the disk. Also, a bar is 
expected to destroy any metallicity gradient within a certain radius, as is observed (e.g. 
G06). Harris & Zaritsky (2009) recently found that the chemical evolution history for 
the LMC is qualitatively similar to the one derived for the SMC (see Harris & Zaritsky 
2004). In both galaxies, the mean metallicity remained at an almost constant, low value, 
from ~10 Gyr until about 4 Gyr ago when the metallicity began a steady ascent to the 
present-day value. In both galaxies, the AMR is consistent with the bursting enrichment 
model proposed by Pagel & Tautvaisiene (1999) (P09). 

7.3. Metallicity Gradient 

The galactic distribution of elemental abundances is central to study the chemical 
evolution of the galaxy. The relative abundances of the elements are sensitive to the star 
formation history, the number of massive stars, the ratio of SN la to II, the importance 
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of AGB stars, the relative yield of the elements, and the exchange of matter between 
the galactic disk, bar and halo through infall or ejection. Measurements of elemental 
abundances throughout the galaxy provide vital input to model the formation and 
evolution of the galaxy. From observations of clusters, 091 and Santos Jr. et al. (1999) 
found no evidence for the presence of a radial metallicity gradient in the LMC. The only 
evidence for a radial metallicity gradient in the LMC cluster system was reported by 
Kontizas, Kontizas, & Michalitsiano (1993) based on six outer LMC clusters (> 8 kpc). 
Hill et al. (1995) were the first to report evidences of a gradient in the field population 
from high-resolution spectroscopy, using a sample of nine stars in the bar and in the disk. 
They found that the bar is on average 0.3 dex more metal-rich than the disk population, 
but the disk stars studied are located within a radius of 2°. Subsequent work by Cioni 
& Habing (2003) also implies a decrease in metallicity while moving away from the bar 
within radius of 6°. Finally, an outward radial gradient of decreasing metallicity was 
also found by Alves (2004) from infrared CMDs using the 2MASS survey. Recently 
Carrera et al. (2008) concluded that the lower mean metallicity in the outermost field 
is a consequence of it having a lower fraction of young stars, which are more metal-rich. 
They found that there is a change in the age composition of the disk population beyond 
a certain radius, while the chemical enrichment history seems to be shared by all fields. 

It is necessary to further investigate the presence, or absence, of an abundance 
gradient in the LMC. We have used the current sample of 101 clusters in this regard. 
Positions on the sky for each cluster are shown in Figure [TH along with the metallicity 
bin into which each cluster falls represented by the color of the plotting symbol. We 
have used the same color scheme as used by G06. Figure [T4T a) represents clusters 
with [Fe/H] determined using CaT method. Star and dot symbols represent data from 
the present work and from G06 respectively. Figure [T4Tb) represents the sample of all 
101 clusters. The adopted center of the LMC (a = 5 h 27 m 36 s ,5 = -69°52'12" ; van 
der Marel et al. 2002) is marked by the cross and the dashed oval represents the 2° 
near-infrared isopeth from van der Marel (2001), which roughly outlines the location of 
the LMC bar. Conversion from right ascension and declination to Cartesian coordinates 
was performed using a zenithal equidistant projection (e.g., van der Marel & Cioni 2001, 
their eqs. [l]-[4]). 

In Figures [15] and [16] we further explore the metallicity-position relationship for 
LMC clusters by plotting metallicity as a function of de-projected position angle and 
radial distance (in kpc), respectively. The symbols are the same as in Figure [[3] 
We can very easily see the tight distribution of the clusters with CaT-based [Fe/H] 
determination. This is an important factor in constraining the SFH of the galaxy. We 
have corrected for projection effects by adopting 34°. 7 as the inclination and 122°. 5 for 
the position angle of the line of nodes of the LMC (van der Marel & Cioni 2001 and 
G06). In this rotated coordinate system, a cluster with a position angle of zero lies along 
the line of nodes, and angles increase counterclockwise (see G06); Radial distances were 
converted from angular separation to kpc by assuming an LMC distance of (m — Mo) 
= 18.5 (~50 kpc); at this distance, 1° is ~ 870 pc. In Figure [TBI we have over-plotted 
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both the MW open cluster metallicity gradient from Friel et al. (2002; dashed line) and 
the M33 gradient from Tiede et al. (2004; solid line). Neither of these disk abundance 
gradients resembles what we see among the LMC clusters, which shows no evidence for 
a gradient. 

Combined, these three figures illustrate that, similar to what was found by 091 (and 
Geisler et al. 2003, G06), there is no [Fe/H] gradient in terms of either position angle or 
radial distance for the higher metallicity clusters in our sample. It is well known that a 
number of metal-poor clusters ([Fe/H] < —1.5) exist in the inner portions of the LMC 
(e.g., 091), suggesting that neither the metal-rich nor the metal-poor clusters exhibit a 
metallicity gradient. In the case of the LMC, the presence of a strong bar component 
may have diluted the metallicity gradient originally present in the star clusters, leading 
to a cluster population that is well mixed (see G06). 

Geisler et al. (2009) also found that the cluster distribution is similar to what 
has been found for red giant stars in the bar, which indicates that the bar and the 
intermediate age clusters have similar star formation histories. This is in good agreement 
with recent theoretical models that suggest the bar and intermediate age clusters formed 
as a result of a close encounter with the SMC. Our findings also confirm previous results 
which show that the LMC lacks the metallicity gradient typically seen in non-barred 
spiral galaxies, suggesting that the bar is driving the mixing of stellar populations in 
the LMC. 

7.4- Kinematics 

The kinematical properties of the LMC provide important clues to its structure. They 
have been obtained from many tracers such as HI gas (e.g., Rohlfs et al. 1984; Luks 
& Rohlfs 1992; Kim et al. 1998), star clusters (Freeman, Illingworth & Oemler 1983; 
Schommer et al. 1992), planetary nebulae (Meatheringham et al. 1988), HII regions and 
supergiants (Feitzinger, Schmidt-Kaler & Isserstedt 1977), and carbon stars (Kunkel et 
al. 1997; Graff et al. 2000; Alves & Nelson 2000; van der Marel et al. 2002). All 
these studies imply that the LMC is kinematically cold, and must therefore to first 
approximation be a disk system (van der Marel 2004). G06 also found their cluster 
sample has motions consistent with that of a single rotating disk system, with no 
indication of halo kinematics. 

To study the kinematics of the LMC, we followed the approach of G06 and 
Schommer et al. (1992) using star clusters as a tracers. To characterize the rotation of 
their clusters, Schommer et al. (1992) fit an equation of the form 

V(6) = ±y m {[tan(0 - O ) sec if + l}" - 5 + V sys 
to their radial velocity data using a least-squares technique to derive the systemic 
velocity (V sys ), the amplitude of the rotation velocity (V m ), and the orientation of 
the line of nodes; they adopted an inclination of 27°. Their best-fit parameters give 
a rotation amplitude and dispersion consistent with the LMC clusters having disk like 
kinematics, with no indications of the existence of a pressure-supported halo. 
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In Figure [T7] we have plotted galactocentric radial velocity versus position angle on 
the sky for our sample (8 from present work + 28 from G06), along with velocity data 
for all clusters listed in Schommer et al. (1992). To be consistent with the approach of 
Schommer et al. (1992), we have adopted the galactocentric velocity corrections given by 
Feitzinger & Weiss (1979). In addition, for this Figure only, we have adopted their LMC 
center (a = 5 h 20' m 40 s , 5 = — 69°14'10" [J2000]) for use in calculating the position angles 
of our clusters (see also G06). We have used the standard astronomical convention in 
which north has a position angle of zero and angles increase to the east (G06); Data 
from Schommer et al. (1992) are plotted as open circles, present data are plotted as star 
symbols and data from G06 is given by solid points. Over-plotted on this figure (dashed 
line) is the rotation curve solution number 3 from Schommer et al. (1992). Derived 
velocities for the present eight target clusters show that their motions are consistent 
with the findings of Schommer et al. (1992) in that the LMC cluster system exhibits 
disk-like kinematics that are very similar to the HI disk and has no obvious signature of 
a stellar halo. Carrera (2008) argued that failure to find evidence of a hot stellar halo is 
related to a low contrast of the halo population with respect to that of the disk, even at 
the large galactocentric radius of their outermost field. Curiously, both old, metal-poor 
clusters as well as intermediate age, metal-rich clusters share the same disk rotation. 

8. Summary and Conclusion 

In this paper, we have studied the properties of CMDs, age, metallicity and radial 
velocities of eight LMC massive clusters. Several of these clusters have no modern 
metallicity determination; out of eight, four have a spectroscopically determined 
metallicity value given in 091 with a larger error. We have used on-average 10 cluster 
members to determine the mean metallicity which shows a scatter of only 0.05 - 0.12 
dex. To derive the metallicity, we obtained medium-resolution spectra at the near- 
infrared CaT using the FORS2 multiobject spectrograph at the 8.2-meter VLT/UT1 
with high signal-to- noise (20 < S/N < 100) and then used the empirically developed, 
simply calibrated method which can make an efficient estimate of metallicity ([Fe/H]) 
for individual RGB stars using the strength of the CaT lines. We also studied in detail 
the CMDs of all eight clusters based on photometric data taken from FORS2 + VLT 
under excellent atmospheric conditions (seeing < 0.7 arcsec) and determined their ages, 
making use of existing deep HST data for several of the more crowded fields. The strong 
CaT lines are also used to determine the radial velocity of each target. 

We have used the derived age, metallicity and radial velocities of all eight clusters 
to study the age metallicity distribution, their relationship and significance, their 
distribution in the LMC and the kinematics of the LMC. We have confirmed the 
efficiency and efficacy of the CaT method in deriving [Fe/H] with much less scatter 
than previous studies, yielding more conclusive results. These data can further be 
used to calibrate the RGB slope-metallicity-age relation based on our homogeneous 
spectroscopic and infrared data. The main results of this paper are as follows: 
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• We have studied the CMDs and reported the age of all eight clusters. Six have ages 
between 0.8-2.2 Gyr and the other two, ie. NGC 1754 and NGC 1786, are very old. 

• We have derived a spectroscopically determined metallicity for all eight clusters, 
with comparatively less scatter than previous studies using the CaT method. Out 
of eight, four have no previously spectroscopically determined metallicity value and 
the rest have large errors in their previously determined values. The cluster Hodge 
2 has no previously reported metallicity value. 

• The metallicity of the six intermediate age clusters, having a mean age 1.5 Gyr, 
is —0.49 with a scatter of 0.04, in excellent agreement with the previous finding 
of G06, which shows a very tight distribution with mean metallicity of —0.49 
(a = 0.09). Similarly, Mucciarelli et al. (2008) used high resolution spectra of 
the giants in 4 LMC globular clusters and found a mean [Fe/H] = —0.38 ± 0.09. 
This agrees well with the theoretical work of Bekki et al. (2004), which indicates 
that a close encounter between the LMC and SMC caused not only the restart of 
cluster formation in the LMC but the generation of the central bar (mean [Fe/H] 
= -0.37 ±0.15; C05). 

• We find that there is no evidence of a metallicity gradient in the LMC, in contrast 
with the stellar populations of the MW and M33, which show that the metallicity 
decreases as galactocentric distance increases. The LMC's stellar Bar may be the 
responsible factor for the dilution of any kind of gradient in the LMC. This result 
is in good agreement with previous work done by G06. 

• We find that our derived cluster velocities are in good agreement with the results of 
Schommer et al. (1992) and G06, confirming that the LMC cluster system exhibits 
disk-like rotation with no clusters appearing to have halo kinematics. 
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Table 2. Position and measured values for all RGB stars. 
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1 
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7.6 


-1.65 


62 


5.24 
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04 
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7.7 
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46 


7.15 
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.24 
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41 
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— 
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1 
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-2.24 
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6.41 





.25 


— 


.65 
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NGC1754-10 


04 


54: 12.80 
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:26:54. 


2 
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7.5 


-1.35 
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6.96 
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.15 


— 
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0.11 
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NGC1754-20 


04 
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3 


38.3 


7.5 


-1.92 
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5.40 
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27 
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04 


54:20.83 
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:25:45. 


8 
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-1.45 


57 


6.14 
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3 
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— 
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9 
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-3.00 


52 


7.62 


5.43 
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.17. 


41 
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7 
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7.6 


— 2 


.27 
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7.47 





12 
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26 


0- 


12 


n 
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05: 


.17: 


:36 


04 


— 69 


:39: 
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5 
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7 


7.5 


_2 


08 


56 
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7.57 


o 


14 


— 0. 


22 


0. 
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05: 


:17: 


■38 


40 


— 69 
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3 


8.6 
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17 


54 
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0. 


11 


n 
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05 


.17 


51 


68 
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7 
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7.6 
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27 


0. 
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05: 


:17: 


.41 


.80 


— 69 
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4 
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3 


7.6 


— 


64 


38 
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7.53 





.13 


— 0. 


24 


0. 


12 


n 


Hodgc02-08 


05: 


.17 


■41 


.62 


— 69 


:39: 
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7 
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2 


10.7 
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42 
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6.53 





.06 
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.60 


0. 


10 


n 
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05 


:17: 
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.16 


— 69 
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3 
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8 


8.9 


_ 1 
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.13 
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24 


0. 
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n 
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05 
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44 


53 
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43 


4 
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9 


7.6 
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6.78 





.01 
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n 
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05 


17 
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37 


7 
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g 
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— 
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7.14 
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07 
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63 


0. 
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05 


17 
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50 
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:38: 
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2 
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8 
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_1 
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43 


8.75 
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.16 
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19 


0. 


13 


n 
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05 


17 


49 


96 


— 69 
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Q 
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9 


7.9 


_2 


.05 
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6.39 


4.90 
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_1 


19 


0. 


15 


n 
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05: 
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49 
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o 
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7.6 


_1 
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0. 
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51 


.89 
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_1 
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.23 
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94 


0. 


12 


n 
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.17. 


■ 5 6 
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2 
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3 


7.5 
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o 


30 
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o. 
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n 
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05 
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36 
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— 69 
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4 
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.14 


40 


8.45 


7.62 





15 
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21 


0. 


12 


n 
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05: 


.17. 


;43 
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— 69 
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3 
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1 


7.5 


_1 


.49 


57 


6.08 


4.99 
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16 


0. 


15 


n 


Hodgc02-27 


05 
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— 69 
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1 
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7.5 
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47 


8.16 


7.34 





.10 
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31 


0. 
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n 
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05 


17 


40 
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08 


4 
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9 


7.6 


_1 


.94 
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7.79 





.18 
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15 


0. 


13 
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05 


.17 
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64 


— 69 
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57. 


9 
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9 


7.5 


— 
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42 


7.11 


6.41 
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— 0. 
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n 
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05 
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.27 


05 


— 69 
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28 


1 
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2 
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— 
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35 
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— 0. 


50 


0. 
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05 
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29 
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— 69 
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2 
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2 
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43 
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o 
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0. 
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05 
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oo 


29 
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g 
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54 


0. 


10 


m 
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19 


14 
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— 63 
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7 
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6 


8.1 


— 


21 


42 


6.49 
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67 


0. 


09 


m 
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09 


13 


— 63 
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5 
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41 


57 
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6.53 


o 


04 
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0. 
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m 
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19 
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56 
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38. 


3 


309. 


2 
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7.67 


6.93 


o 
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0. 


11 
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05 


19 


08 


82 


— 63 
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3 
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0. 
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05: 


19 
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97 


— 63 
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50. 
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8 
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57 


6.72 
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.03 


— 0. 
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o. 


09 


m 


NGC1900-21 


05 


19 
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35 


6 
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2 
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— 
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36 


0. 


13 


m 
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04 
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17 


6 


296- 
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54 
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6.46 
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0. 


10 


m 


NGC1900-26 


05 
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9 
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7.74 
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.09 
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53 


0. 


11 
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19 


09 


.75 


— 62 
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2 
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3 
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6.86 
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— 0. 
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0. 


11 


m 
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19 


17 


.94 


— 63 
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1 


333 


3 
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46 
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39 


o, 
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n 
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05 
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11 
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58 


0. 


10 


n 


NGC1900-09 


05 
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.4 
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9 
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6.90 
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09 
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30 


4 
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39 
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— 0. 
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0. 
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19 


06 
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o 
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05 
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29 
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9.0 


_1 


07 
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.30 


_1 


15 


0. 


13 


n 
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42 
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1 
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6.01 


4.94 
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_1 
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05 
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_1 
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_1 
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0. 
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05: 
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g 
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40 
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29 


46 
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47 


2 
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0. 
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05 


29 


46 
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5 
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6.75 
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52 
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1 
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53 
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6 
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7.6 


— 
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6.72 
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0. 


10 


m 
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29 
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Table 2 (cont'd) 
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Star ID Q2000 ^2000 KV tr RV V — Vh B S/N W W 1 a w , [Fe/H] ^[Fe/H] Comment 

(h:m:s) (d:m:s) (km s _1 ) (km s _1 ) (mag) (dex) (dex) 
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Table 3 Derived LMC cluster properties. 



Cluster 


Age 


N 


RV ± ctrv 


RV (091) 


[Pc/H] ± <7 [Fo/H ] 


[Fe/H](091) 


[Fc/H] (091 6 ) 


[Fe/H] ± ^fFc/Hl (photometric) 




(Gyr) 




(km s _1 ) 


(km s _1 ) 


(dex) 


(dex) 


(dex) 


(dex) 


Reference 


NGC1795 


1.3 


9 


255.0 ± 6.8 


266 


-0.47 ± 0.10 


-0.23 


-0.33 






NGC 1754 


14.0 


8 


234.1 ± 5.4 


236 


-1.48 ± 0.09 


-1.54 


-1.28 


-1.42 ± 0.15 


Olsen et al. (1998) 


NGC1786 


12.3 


7 


279.9 ± 4.9 


264 


-1.77 ± 0.08 


-1.87 (-1.75 ± 0.01 a ) 


-1.59 


-2.10 ± 0.3 


Brocato et al. (1996) 


Hodge 2 


1.6 


7 


266.3 ± 9.7 




-0.49 ± 0.12 










NGC 1900 


0.8 


10 


305.0 ± 8.3 


300 


-0.55 ± 0.10 


-0.81 


-0.70 






SL509 


1.2 


13 


317.0 ± 8.4 




-0.54 ± 0.09 






-0.65 ± 0.2 


Piatti et al. (2003) 


SL817 


2.2 


10 


281.1 ± 4.5 




-0.41 ± 0.05 






-0.35 ± 0.2 


Piatti et al. (2003) 


SL862 


1.7 


13 


238.0 ± 6.2 




-0.47 ± 0.07 






-0.75 ± 0.2 


Piatti et al. (2003) 



a: Based on high resolution spectroscopy: Mucciarclli et al. (2009) 
b: Transformed using equation in text 



Table 4 Metallicities of young, intermediate and old clusters. 



Age range 


CaT data 




Literature data 




all 


data 




(Gyr) 


[Fe/H] (dex) Agc(Gyr) 


N 


[Fc/H] (dex) Agc(Gyr) 


N 


[Fe/H] (dex) 


Age(Gyr) 


N 


< 1 


-0.55 ± - 0.8 ± - 


1 


-0.34 ±0.18 0.2 ±0.2 


35 


-0.35 ± 0.18 


0.2 ± 0.3 


36 


1-3 


-0.47 ± 0.09 2.0 ±0.5 


26 


-0.64 ±0.18 1.8 ±0.4 


22 


-0.55 ± 0.16 


1.9 ± 0.5 


48 


> 12 


-1.65 ± 0.22 14.1 ± 1.3 


7 


-1.58 ±0.29 14.5 ± 1.4 


6 


-1.62 ± 0.26 


14.3 ± 1.3 


13 
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Figure 1 Sample spectra from RGB stars in our target clusters covering a range in 
metallicities. The three Ca II lines are marked for reference; the change in CaT line 
strength with [Fe/H] is readily visible. Calculated summed equivalent widths and 
metallicities for each star are given. 



28 








f-J --4-i- 




i 





300 600 900 1200 1500 

X position (pixels) 



100 200 300 400 500 600 

Distance (pixel) 







350 










300 










250 


T 










/ zfe \ 

/ ♦ 


trie Velocity (km/s 












\. # V 


Heliocen 















-50 




I 







900 

X position (pixels) 



300 400 
Distance (pixel) 



Figure 2 (left panels) The (x, y) positions for the stars in cluster NGC 1754 (upper) 
and NGC 1900 (lower). The adopted cluster radius is marked by the large open circle, 
and stars outside this radius are considered non-members. The blue color triangles 
are the stars which are selected from the radial velocity criterion and then the final 
cluster members (purple inverse triangle) from these blue stars have been selected on 
the basis of metallicities (see text for details), (right panels) Radial velocities for our 
spectroscopic targets as a function of distance from the NGC 1754 (upper) and NGC 
1900 (lower) cluster centers. The horizontal lines represent our velocity cut and have a 
width of ±11 km s _1 . The cluster radius is shown by the vertical line. The error bars 
represent the random error in determining the radial velocity for each star. 
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Figure 3 Metallicities as a function of distance for all target stars in NGC 1754 (left) 
and NGC 1900 (right). The [Fe/H] cut of ±0.20 dex is denoted by the horizontal lines. 
For NGC 1754, a metal-poor cluster, the field is relatively easily distinguished from 
the cluster. The plotted error bars represent the systematic error in [Fe/H] calculation. 
Triangles are stars satisfying the radial velocity criteria and the star symbols are those 
which qualify both selection criteria based on radial velocity and [Fe/H]. 



12 - 
10 - 




2 - 
- 

0.0 -0.5 -1.0 -1.5 -2.0 -2.5 

v-v HB 



Figure 4 Summed equivalent width vs. brightness above the HB for all stars considered 
to be members of NGC 1900 (circle) and NGC 1754 (triangle). The slanted lines are 
iso-abundance line with a slope (3 = 0.73 at the mean metallicity of the clusters i.e. 
[Fe/H] = -0.55 and -1.48 for NGC 1900 and NGC 1754 respectively. 
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Figure 5 CMD for the entire NGC 1754 (left) and NGC 1900 (right) field, with CaT 
stars marked with big dots. Confirmed cluster member RGB stars are also shown by 
star symbols. 




Figure 6 CMDs of NGC1795; NGC1754; NGC 1786; Hodge 2; NGC 1900; SL 509; 
SL817; and SL862. The grey dots represent all stars from our photometric lists, the 
most probable cluster members (see the text) are plotted with dark dots. The dashed 
and solid lines show the MACHO mean RGB region and the location of the field RR 
Lyrae star, the big square shows MACHO RC region. 
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Figure 7 The stellar surface-density distributions (number of stars per arcmin 2 ) of 
NGC1795. Left panel: All stars in our photometric list. Right panel: Most probable 
cluster members. 




Figure 8 The statistically decontaminated CMDs of NGC 1795; NGC 1754; NGC 1786; 
Hodge 2; NGC 1900; SL509; SL817; and SL862. The dark triangles represent 
spectroscopically observed stars. The solid lines show best fit of the Padova set of 
theoretical isochrones (Marigo et al. 2008), whose age and metallicity are marked. The 
smaller panel for NGC 1754 represents the HST VI data for the cluster. 
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Figure 9 Left panel: CMD of SL 509. The regions of the two suspected RCs are marked 
with large squares. The solid lines show the Padova theoretical isochrones (Marigo et al. 
2008) for ages of 1.0, 1.2 and 1.4 Gyr. Middle panel: The luminosity function of the RC 
region. Right panel: The positions of the main RC stars (filled circles) and suspected 
secondary clump (crosses). 
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Figure 10 Difference between present and previously determined [Fe/H] value (A Fe/H = 
Fe/H present - Fe/Hnterature) vs. Age diagram. Filled and open circles represent original 
and transformed 091 values respectively and open triangles are for photometrically 
observed clusters. Cross marks the high resolution value given for NGC 1786 by 
Mucciarelli et al. (2009). 
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Figure 11 Age distribution of LMC clusters. 
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Figure 12 Metallicity distribution of LMC clusters. 




Figure 13 The Age-Metallicity Relation for the LMC. The points with error bars are 
101 LMC star clusters which have age and metallicity measurements in the literature 
and in the present study. Star symbols are the points from the present study, closed 
circles from G06 and open circles from Harris & Zaritsky (2009), respectively. Pagel & 
Tautvaisiene (1999) burst model is shown as a solid line. 
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Figure 14 Positions on the sky and derived metallicities for our target clusters. 
Metallicity bins are given in the lower left corner of the plot. The adopted LMC center is 
marked with the cross, and the dashed line roughly outlines the bar. Panel (a) represents 
clusters with [Fe/H] determined using CaT method, symbols are as same as in Figure 
[TBI Panel (b) represents the sample of all the 101 clusters. 
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Figure 15 Metallicities of our target clusters as a function of de-projected position 
angle. We have used the LMC geometry of van der Marel & Cioni (2001) to correct 
for projection effects. This plot illustrates that there is no apparent relation between 
position angle and metallicity in the LMC. Symbols are as same as in Figure | 
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Figure 16 Cluster metallicities as a function of de-projected distance (in kpc from the 
center of the LMC). We have assumed a distance of (m — Mo) = 18.5. Over-plotted 
are the metallicity gradients observed in the MW open clusters (dashed line; Friel et 
al. 2002) and M33 (solid line; Tiede et al. 2004), which help to further illustrate that 
the LMCs cluster system lacks the metallicity gradient typically seen in spiral galaxies. 
The lack of a gradient is likely caused by the presence of the central bar (Zaritsky et al. 
1994). Symbols are as same as in Figure [TBI 
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Figure 17 Galactocentric radial velocities as a function of position angle on the sky 
for the clusters in our sample (stars), from G06 (filled circles) as well as those from 
Schommer et al. (1992; open circles). Rotation curve solution number 3 from Schommer 
et al. (1992) is over-plotted as the dashed line, showing that all data sets are consistent 
with circular rotation. 



